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PREFACE 

The Agriculture and Resources Inventory Surveys Through Aerospace Remote 
Sensing Is a multiyear program of research, development, evaluation, and appli- 
cation of aerospace remote sensing for agricultural resources, which began In 
fiscal year 1980. This program Is a cooperative effort of the U.S. Department 
of Agriculture, the National Aeronautics and Space Administration, the National 
Oceanic and Atmospheric Administration (U.S. Department of C> nerce), the 
Agency for International Development (U.S. Department of State), and the 
U.S. Department of the Interior. 
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1. INTRODUCTION 


Multi frequency microwave radiometer (MFMR) measurements were taken by the 
National Aeronautics and Space Administration (NASA) C-130 aircraft at 1000 
and 1500 feet above ground level as part of the Agricultural Soil Moisture 
Experiment (ASME). This experiment was conducted near Colby, Kansas, on 
July 18, 20, 21, and 22, 1978, and on August 8, 9, and 11, 1978. These data 
are a measure of the natural microwave emission of the targets as seen by the 
sensor. The MFMR Is a group of sensors that collect data at five frequencies: 
1.42, 5.0, 18.0, 22.5, and 37.0 GHz. Location of these sensors on the 
aircraft Is shown In figure 1-1. Each sensor collected data by viewing the 
ground with either a 0° or 40° Incidence angle; only one Incidence angle could 
be used at any one time. Because of mechanical Incompatibility of their 
antennas, the 1.42 and 5.0 GHz frequencies could not be used simultaneously. 

Soil samples from several layers were collected from preselected fields of 
approximately 40 acres on each of the 7 days of flight. These samples were 
weighed, oven-dried, and weighed again so that the moisture content of the 
layers could be calculated. The moisture contents will be used for comparison 
with the MFMR data. 

Simultaneously with the MFMR data collection, the Precision Radiation 
Thermometer Model 5 (PRT-5) measured the Infrared temperature of the emitting 
surface. Knowledge of this temperature may allow the conversion of the MFMR 
antenna temperatures to emlsslvltles. 

A method of converting the and PRT-5 computer compatible tape (CCT) data 
Into disk files Is outlined In this document. Each disk file will contain the 
date, the MFMR data, the PRT-5 temperature data, and the ground reference posi- 
tion of the data within the sampled field. This conversion Is accomplished by 
four processing programs. These ore listed and each program Is discussed In 
this report. 
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Flgurt 1-1.- Botton view of the NASA C-130 aircraft. 



2. REQUIRED INPUT DATA 


As the aircraft flew down the flight line, the WMR and PRT-5 data were 
recorded on tape In an analog format. The tapes were sent to the Sensor 
Analysis Laboratory (SAL) located at the Lyndon B. Johnson Space Center (JSC), 
Building IS. At the SAL, the analog data were digitized and recorded on tape 
along with the Auxiliary Data Annotation System (ADAS) data time of acquisi- 
tion. Each record In the file also contains the roll, pitch, and drift angles 
along with the altitude, ground speed, and true heading of the aircraft. The 
data and flight parameters were measured approximately every 0.6 second. 

A great deal of aerial photography was collected as part of the ASME. 
Photographs taken at an attitude of 8000 feet were used to construct con- 
trolled strip mosaics of each flight line. Additional aerial photography was 
acquired as the aircraft collected the MFMR data. The acquisition time and 
frame number of every photograph were recorded on the analog data system. 

This made It possible to determine the aircraft's position at the frame times. 
The camera positions and frame numbers were plotted on transparent overlays by 
the JSC Cartographic Technology Laboratory. Additional overlays were made 
showing the location of the sampled fields In each flight line. All overlays 
were made at the same scale as the controlled strip mosaics. 
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3. PROGRAM EXPLANATION 


3.1 CONVERT PROGRAM 

The MFMR data as provided by the SAL are grouped Into records with a length of 
160 characters. Some characters In each record are written In the American 
Standard Character Code for Information Interchange (ASCII) character code. The 
rest of the characters are In binary form. The computer used In this study Is a 
National Advanced Systems AS/3000 which uses the Extended Blnary-Coded-Declmal 
Interchange Code (EBCOIC) for character representation. Therefore, before the 
ground registration of the data can begin, the ASCII and binary characters must 
be converted to EBCDIC. This Is accomplished by a FORTRAN program called 
CONVERT. A listing of the program and Its execute file are given In appendix A. 
The program reads In the data, does the character conversion, and reformats the 
data to a record length of 132 characters so that a printed copy can be 
obtained. The data file from the SAL can contain three types of records: call* 
bratlon, baseline, and data. The program CONVERT will write each to a separate 
file In EBCDIC representation. A sample listing of the output data file Is 
shown In appendix A. The column titles are added for clarity but are not 
actually written to the data file. Note that at this point all columns have 
numbers In them, ev*n If the sensor was turned off during the run. To Identify 
the useful data, che file Identification code must be referenced. Only the 
useful data will be placed In the file output by the final processing program. 

3.2 MFMR PROGRAM 

The second program used to process the MFtffi and PRT*5 data Is a FORTRAN program 
called MFMR. A listing of the program and Its execute file are given In 
appendix A. This program reads In the reformatted and converted data file and 
computes the location of the aircraft's negative z-axls Intersection with the 
ground In a scene-based coordinate system. The Iff MR program refers to the 
camera location but, as clearly shown In figure 1-1 of section 1, the MFMR 
sensors have an along-track displacement from the camera. The sensor displace- 
ments and beamwldths are listed In table 3-1. To obtain the location of the 
sensor footprint center, the displacement must be added to the along-track 
displacement. 
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TABLE 3-1.- MFMR SENSOR SUMMARY 


Sensor 

Along- track 
displacement 

(ft) 

Beamwldth, degrees 

Band 

Freouency 

(GHI) 

Wavelength 

(cm) 

Half- power 

Null to null 

L 

1.4 

21.2 

+25.0 

16.0 

40.0 

C 

5.0 

6.0 

+25.0 

4.5 

12.0 

Ku 

18.0 

1.66 

-42.0 

5.0 

14.0 

K 

22.0 

1.36 

-42.0 

5.0 

14.0 

Ka 

37.0 

0.81 

-42.0 

4.0 

12.0 


Three Inputs from the terminal are requested by the program. The first Input, 
called AMISS, Is the distance (In feet) of the aircraft toward the north from 
the southern field boundaries If the flight line runs east-west. If the 
flight line runs north-south, then AMISS 1$ the distance toward the west from 
the eastern field boundaries. AMISS Is measured at the beginning of the 
flight line. The second Input, called YUP, Is the crosstrack distance (In 
feet) that the aircraft's position changed between the beginning and the end 
of the flight line, rlgure 3-1 Is a diagram of a flight line and shows the 
distances represented by AMISS and YUP. Both Inputs are measured with the 
overlay on the strip mosaic. The final requested Input Is a three-symbol 
numeric Identifier called Code. The Codes used for the ASME MFMR data are 
shown In table 3-2. 

TABLE 3-2.- IDENTIFICATION COOES FOR THE COLBY ASME DATA 


Code 

Day, 1978 

Sensor 

Polarization 

1 

199 (July 18) 

13.3 GHz Scatterometer 

Like 

2 

201 (July 20) 

4.75 GHz Scatterometer 

Cross 

3 

202 (July 21) 

1.6 GHz Sc attend* +.er 

Horizontal receive 

4 

2C3 (July 22) 

0.4 GHz Scatterometer 

Vertical receive 

5 

220 (August 8) 

0* L, Ku, K, Ka MFMR 

Horizontal and vertical receive 

6 

221 (August 9) 

40* L, Ku, K, Ka MFMR 


7 

223 (August 11) 

0* C MFMR 


8 


40* C MFMR 
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Program MFMR creates an output file that contains the downtrack and crosstrack 
locations of the airplane with respect to the beginning of the flight line. 
This file also contains the corresponding MFMR data values for those locations 
and the time that the data were acquired. Time of data acquisition Is no lon- 
ger needed for MFMR data ground registration, but It will be used to register 
the PRT-5 data acquired simultaneously. 

3.3 MPLOT PROGRAM 

The third program, W>L0T, Is a Statistical Analysis System (SAS) program which 
reads the file output by program WMR and plots the aircraft's ground track 
along with one band of corresponding MWt data. The plot Is at the same scale 
as the strip mosaics. A program listing Is given In appendix A, and a portion 
of the plot is shown in figure 3-2. At this point In the analysis, the MFMR 
data are referenced by the distance (In feet) downtrack and the distance from 
the southern or eastern field boundaries, depending upon the direction of the 
flight line. It Is necessary to know the ground reference position of the 
data within the sampled fields. This Is accomplished by using the overlays in 
conjunction with the plot. Both overlays, one with photographic position and 
the other with field boundaries, are placed on the plot in the following 
manner. First, a time is found when the aircraft's flight parameters are 
available and when an aerial photograph was taken. The time represented by 
each asterisk In the flight -path plot can be found by using the exclamation 
points plotted alongside. The exclamation points are time marks when the 
aircraft clock was at the minute or a multiple of 10 seconds after the minute. 
Next, the overlays are placed on the plot so that the photographic position 
from the overlay Is on top of the asterisk representing the same time. The 
solid line paralleling the plot of the flight path should be even with the 
southern field boundaries for an east-west flight or with the eastern field 
boundaries for a north-south flight. Figure 3-3 is an Illustration of the 
plot with the overlays In place. The downtrack distance (In feet) of each 
field's closest boundary to the beginning of the flight line Is read from the 
plot. These distances, along with the dimensions of the corresponding fields, 
are written In a separate file. This file Is used as an Input to the next 
program. 
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Figure 3-2.- A portion of the plot of the MPLOT program 




When the procedure Is carried out In the manner described above, a dl crepancy 
may become apparent. If the overlay and plot were registered at an ext erne 
end of the flight line, the overlay and plot positions may not match at the 
opposite end of the line. This Is caused by a lack of sufficient accuracy In 
the recording of the aircraft's flight parameters. Therefore, It Is recom- 
mended that the overlay and plot be registered near each sampled field before 
the downtrack distances are read. 

3.4 MGRIO PROGRAM 

The final processing program, called MGRIO, reads In the boundary data file, 
the file created by program MFMR, and the corresponding file of PRT-5 data. A 
listing of program MGRIO and Its execute file are given In appendix A. 

Program MGRIO determines which MRtt data lie within a sampled field and 
calculates the location of the sensor footprints with respect to the northern 
and western field boundaries. The program then searches through the PRT data 
to find the infrared temperatures that correspond to the footprints. One 
output file Is created for each sampled field within the flight line. The 
output files can then be combined In the manner which best suits the analysis 
technique that will be used. 
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PROGRAM AND DATA LISTINGS 


APPENDIX A 


PROGRAM AND DATA LISTINGS 

The following data are presented in this appendix. 

a. CONVERT FORTRAN 

b. CONVERT EXEC 

c. Converted and reformatted data listing (See figure A-l. 

d. MFMR FORTRAN 

e. MFMR EXEC 

f. MPLOT SAS 

g. MGRID FORTRAN 

h. MGRID EXEC 
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FILE I CONVERT FORTRAN A EOOL / JOHNSON SPACE CENTER 
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THIS program 
II NARY AND A 

three eiles 

data, caubra i I un »nu u*btl_iric. 

WRITTEN BY 1 JOHN C RICHTER FEBURARY. 1981 


READS JN ONE DISK FILE OF MFMR DATA, THE 
BINARY AND ASCII CHARACTER COOES ARE CONVERTED TO EBCDIC. 
THREE FILES ARE WRITTEN, ONE FOR EACH COLLECTION MOOE ! 

' RATION AND BASELINE. 


cccccccccccccccccccc 

I OgfCALM A (10 
NTEGER , i B (28 
NTEGER OCHARl 
AT* DCMAP/C4, 


CCCCCCC5S5SCCCfCCCCCgCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
... ... 1*951 .CM* 

_ 5CHAP/ 64,64,64, 64 .64, 64, 6*. 64, 64, 64, 64, 64, 64, 64, 64, 

164 .64.64.84 *64 .64 . 64.64*64 .64 .64.64 .64.64.64 *64.64 . 

290. 127. 123. 91,. 08. 80. 125. 77. 93. 92, 78. 107. 98. 75. 97. 240. 241, 
3242. 243. 244, 245. 246. 247. 248. 249. 122, 94. 76. 126. 110.11 1. 
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FORMAT (108AI.26A2) 
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A(K)4DCHAR)A(K) ) 
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CMX*A( 108) 


101 

10 


t 
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10 K*|.l 

[•DCHAR? 
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214. 13. IX. 14. 1 X. 14. 13. 14. 1 1. IX, 14, 1 3.1*. 13. IX. 14. 1 3. I* 
2 CONTINUE 
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2) .8(3) .0(4) .0(7) .0(8) .8(91.81 10) , 


IX, 14,13, 
4,13) 
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N00 130 
CONOOlAO 
CON00150 
CONOOIOO 
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00090 0 9999009 «■ * « ■ »' ■ - ■ — — - a h * 1 

4i4i4i4i44ir»i4t44i4t/HA4i4t/H/H™/H4^^ , 

OOOOOOOOOOOOOOOOOOOOOOOOOOOO I 

000000440000000000000000044 ® ) 


r^nrwnn«r>nw>*nini*^n«vnnf^nr>fnr>mr>r>«nr» 


CM 
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Figure A-l.- Sample output from program CONVERT and an explanation of column values. 





















ORIGINAL F' ~ r 
OF POOR QL= • 1 -‘ - • 


F!l f : MFMR 


FORTRAN A EOOL / JOHNSON SPACE CENTER 


X.F5.1.1X.F5.1. 


KATMPV(l).KASDVU) . 

format II3.H.|1.I1.*x.I 2.|X.F4,I , 1X.I4.1X.F5. 1*1 
IX.FS.I.lXtFS.lflXf! JjRXjIA.p.lXjlA.M.lA.I j. 
&X.I4.13.I4.I 3. IX. 14.13.14. i3.lX.l4.I3.i4.i3) 

SPEEO I I 1 ‘FLOAT • iSPCCO) 

ALT»FLOAl HALT) 

MlN*FLOAT ( IMIN) 

ACSEC (II* (MIN‘60,0) *SEC 
OT0»« (ROLL* 3.1 A 1543) / 180.0 
YOFF ( I)«AlT‘TAN(OTOR) 

Ni|MBAC*I-l 
I MIN* I MI N* 1 000 

Isec»!f|x isecmo.o) 

ITME<|)«INlN‘fSEC 
SPEED 1 1 ) *SPEEO ( I ) • 1 . 68T8 
■■ hEA0.lt. *n.O)HEAD«HEAO*3«>0.0 


cIr»I-1 

re — 

F ( 

F ( 


.NF.nICHK.ITIMEtll-lTIMFtICTR) 

i CMK.GT.30)IFLAO«l 
CHK.LT. 0) JFLAG-? 

FLAG.EG.1)ITTmE<II‘ITIMF<I)- 4-- 
FLAG.EO.?ilTlME(j)-}TjME(j 
FLRG.E'J't* XCSEC ll)‘*C5EC(I>.3 
ack(Ii«ooIftiI)*heaO 


i-aoo 
►60000 
►3600.0 

1 TftAC" : 

?0) CONTINUE 

JF(ILINE.F'J.1)CH8SN«-1.0 

|f(!lJne.eo.a)chgsn«-|.o 

!f(Il!ne.fo.s)chgsn«-i.o 

YUP*YUP/ ( I ACSFCINUMBAC) -ACSEC 111) ‘SPEED U » ) 
NN‘NUMRAC‘1 

these STATEMENTS COMPUTE the AVERAGE FLIGHT, DIRECTION 


1 A 


ANNsFLOA T (NN) 

00 1* K»1,NN 
hTOT«HT0T*TRACK(K) 
havg«htot/ann 


this SECTION COMPUTES THE LOCATION OF THE AIRCRAFT 


00 11 I.l.NIIMRAC 

ELPTFE*FL0AT(ITIME(I*1)-ITIME(I) 1/10.0 
1F(ELPTME.GT.A0,0IELPTME»ELPTmE-*0.0 
MSEC*IF I X tELPTMt ) 
IF(MSFC.OT.Z)ELPTME=ELPTME-msEC 
angle. I ( (HAV fi-TRACK (I ) ) >‘3. 1*1 S9 31/180.0 
01 SP*ELPTME ‘SPEED < I ) 


x6 

TO 

L* 

YD 

TO 


SP‘01 SP‘COS ( ANGLE ) 

SP.OISP‘STn<ANGLE)‘CHGSN 
♦1 

Ftl.FQ.llGO TO 1Z 
FF«<YOFF(f)-YOFF( 1-1) ) ‘CHGSN 

TOT? (U*TOT? { L - 1 > •?§ISP*Y0IFF» (YUP‘TOTX (L) - YUP‘TOTX (L-l ) ) 
GO TO 11 
i? TOTX (L) «XD1SP 

TOTy(LI‘YOISP*YOFFU)*AMISS 
11 CONTINUE 

THIS SECTION DETERMINES WHICH DATA IS AT 10 SECONDS 


00 


Kid 


8 I.L.NIIMRAC 
1-1 

EN( I ) .0.(1 

F(iniFF.LE.3)ITEN(I)-TOTY(L)*100.0 
0IFF.GT.R6I ITEN( I > • TOT Y ( L > • 100.0 


(EMP‘100) 


DO 


.Ft, _ 

cGnt Inije 
* s K.l.NUMRAC 


MFM00800 

MFM00810 

MFM00820 

MFM00830 

MFM00840 

MFM008S0 

MFM00860 

MFMOORTO 

MFM00880 

MFM00890 

MFH00900 

sssssu 

MFMG0930 

MFM0094Q 

MFM00950 

MFM00960 

MFM009T0 

MFMOO980 


MEMO 

MFMO 


10* 


WRITF(13.10*)!Cl.ICZ.IC3tlLlNE.IRUN,TOTX(K) .TOTYtK) .ITEN(K) , 

LTMP ( k ) .LSO ( K ) .CTMPH (K ) .CSOH (K ) .CTMP vtKi .CSDVtK) .KUTMPH (K ) . 
KUSDHtK ) .KUTMPV(K) .KUSOV (K ) .KTMPH(K) .KSUM(K) .KTMPV (K) ,KSUV (K) 
KaTMPH(K) .KASDH(K) .KATMOV (K) .KASUVtK ) . I T I ME <KI 
FORMAT < - 


STOP 

fnO 


)RMAT ( TI 1 « 1 X . | . •/ • • I 1 . 1X.F8.0. X .f S.O* HfFS.O.ZX. I*. IX. 
M./x.U.lx, 3< IX. I*. IX. 13.ZX. «. IX. 13. jx. 16. IX. I 3.Z> • 
I*. 1*. 13. IX. 4.1X, 13. ?*.!*. IX, 1*1X.I4.1X.13.3X.!6) 


MEMO ISO 0 
MFMO I S 1 0 
.MFM015Z0 
MF MO 1 S3 0 
MFMO 1 S* 0 
MFMOibSO 
MFMO I S60 
MFMO I S TO 
MFM01S60 
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ORIGINAL PACT. i 3 
OF POOR QUALITY 


FILF! MFMR E «FC * EOOL / JOHNSON SPACE CENTER 


HFMOTE TO HOUSTON 
SPOOL F NOHQLO 

C.LOPAL TUTLjS CHSLIB FORTMOO? 
FTLFOEF S RFaDER (PERM 
FILEOEF f, PRINTER (PERM 
FTLFOEF 7 Punch (PERM 

I SK si EOATA A (PERM 
. , - xJS* M 7MFMH « (PERM 

Fjl.fOgF )S TFRMJNAL (PERM 


FTLFOEF I) pi 

n\m ii w 

rfirw »- — 

LOAO MPhS 


iFO^f j* TFWMINAL (PCOM 


STA°T 


g|CFM 


F 

CKH F LPfCL 




HI 


hi 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


nir: nplot s*s * root / Johnson space centeb 


C'*S EICEOEE F ILE 1 DISK HFMB1 YMFMP Al 
OAfA PABTl P*UT?I 
INFILf ElLfll 

n»Uf *0 X 6.0 •IN Y S.o Y T I C S.O «39 YJ 

*A«I lYl/in.OWrO. 01/3,01 
IF YT ic « 0.0 THEN YTIC-.I 
YY»*«»J00. 0*3000. 01 
opop *At 

jF X LE 40000 THFN OUTPUT PAP T 1 I 
IF X r,F ?0000 THFN OUTPUT PAHT7I 
ppor P| or 0 *Ta«PaqT|i 

p| OT X ® Y • * * • X»YT IC«* ! • X • Y • /OVEPl 4Y 
HA»IS« -3000 TO 9000 PY 1000 
v*x: . 0 TO *0000 BY looj 

HSP •’ If 10 VSP4CF* " 
hhFF.- •''00 3150 <4500 
VPFF.O tJ 40000 BY 5000 
VPns* iso I 

P30C PLOT F)*T*«P*MT?I 

P| OT *«Y»>«> X*VTIC«* 1 • X» Y?« •• • /OVEBLAY 
hax)S«-3000 TO O000 BY 1000 
VAX(S*20000 TO N?000 BY 1000 
H5P4CEM 0 VSPACF *B 
HPFF,0 3000 3150 6500 
Yl’tr.^oOOO TO 6?000 BY 5000 
vPOS«350 I 


I 
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ioo^v^nnn<vinnrxvvT,'>nnfv»-i-) mm v 1 


rati MGRIO 


U.v . • 

OF POOR - v u 


F OR T WAN D EOOL / JOHNSON SPACE CENTER 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCl.CC 



ERATUPES APE WRITTEN IN EACH FILE 


the inputs from unit i arei 
ImIN • The number 6f MINUTES PAST THE HOUR ON The 


SEs 


NERDAS CLOCK. 

The NUMBER of SECONDS PAST THE MINUTE ON THE 
NERDAS CLOCK 


PTEMP • THE ARRAY CONTAINING THE PRT-5 TEMPERATURES . 
IRANGE • THE TEMPERATURE RANGE SETTING UN Thf PRT-S 
Tat * The array CONTAINING THE TOTAL AIR TEMPERA 
AS MEASURED AT The AIRCRAFT 


TUBE 


THE INPUT VARIABLES FRON UNIT 9 ABE! 


&P: ?£ ... 

» Th? NLIMHER 

X » THE pOWNTR 
V ■ THE plSTAN 
BOUNDARIES 


HUN 

ME MR 
MF MR 


TABI E IN TEXT FOR EXPLANATION 

flight line numher 

'H OF THE Ruf 

RACK gI|TANCE (IN FEET) OF 


THF DATA 
FIELD 


lSTANCt TO The SouTHERN/EASTtRN 
— ARIES 

IMF MR • THE ARRAY CONTAINING Th£ RADIOMETER UATA 
The INPUTS from UNIT 10 ARC! 

IFLO • AN ARRAY CONTAINING The numbers OF the SAMPLtU 


START 

!S!8E 


EI 5 LDS IN »HE V^HTi.lNE 


TANCES OF THE FIFLOS IN Tnt UNE 

' *” K ) 


ThF OOWNTHACK OlSfANCtS OF THE fiflu 
Thf LENGTH of The fields (alongtralk 
Th£ WIDTH OF The FIELOS (CROSSTRACK) 


The W 

The OUTPUT VARIABLES ARE! 

i C, ’ IC S’i5^S L 8/ l M^*^W^HERS B ^Lo BOUNDARY 
y » NUMBtR OF FEET from Tm£ WESTERN FIELD bOUNUAwt 

MEMR » THE ARRAY OF THE PAD I OMt TER DATA 

the write statement useo is dependent upon the 

DATA Type. MECHANICAL LIMITATIONS PREVENTS i-8»NU ANU 
e-RAND from being ACQUIRED S I MUL T ANEOUSL Y . FOR THE CULbY 
EXPERIMENT. THE KU. K AND KA BANU SENSORS WtRE TuRNtD ON 
WHENEVER The L-BANO WAS ON. The C-BANO sensor was OPtRATEO 
Br ITSELF. 


CCCCCCLCCCC^CCCCCCCCCCCLCLCC 


cccccc S5?raS£S c ??£S??S^fS^5£^S^SS^Miiiis;i«7;*«^;i«i.i*ix ! 6. 

REAL START ( 1 ?) .XWlOE ( IS) .YW|OE ( IS) • APR T (ISI.ATAT(IS) 

NTEG|R HFMRY(^O) «K tJme ( IS) »mFmrx ISIOI • I T Ime (51 0) .PI IMt irooi 


. y (20 


INTEGER 
INTEGER MMFMR 

Common a 

real PTEMP(rOO) .TAT (400) 
MEAD THE INPUT files 


100 


101 


READ (9. 1 001 IC1 . IC2. ICO. 1L INE. IRun.mFMRx ( l ) . mFhry ( 1 ) . 

( ImEmH ( 1 ,N) .N«l. 8 )»ITTm£(1I 
format oi 1 .lx.lj. X.1.X.IS.1X.J5.9X.I-.1X. 

! 3.2x « |4 , j X . 1 J» A,l4.ix,|3.ix.J‘..lX.IJ.Ix,|4.1X.IJ.FX 
4. IX. 1J. lx. 1 *. X, 1J. 2x, 14, lx. lJ.Ix. 14, IX.l3.JX.lt!) 

M«2 «5 1 0 

AO (9*101. 


FORMA 


m .ENO* 
*•15. 


i I ) MF MR X ( M ) • MF MR Y ( M ) 
i*tI5..9X,l4.|X, 


( I MF MR ( M « N ) , N * 1 . IN) . I ! IME (M 


1 I3.?X.l4,lX.l3.1X.14,i*,I3.^X.l4.1*,IJ.l* 

2 i 4 .lx.l 3 . 1 x. 14 .lx. 13.2*, 14. lx 


1 NUMBiM 

)I CONTINUE 


.13. 


. 4.1A.1J.2A, 

X . U , 1 X . J 3. JX . I 61 


DO 2 N»1 .IS 

READ ( 10.102 .END* [2 ) IFLO (Nl .START (N ) .aw lot (N) , »w IDF i 
102 FORMAT (12. lx. FT. 1 . lx. Ft,. l.U, Ft,. 1) 

2 NFLDS»N 
12 00 400 mk k 1 , 4 0 0 

READ (2* 404, t. NO *40 1 ) 1 MIN. s.;c .PTf. HP (MK ] , IRANut . I A I ( MK ) 
404 r OPM AT (?0x, 12, H.F».l,sj«.FS.2.|4«.Al.v».t :,.«!) 


MGROOoio 
mGmOOO 30 

MGR00040 
MGH00050 
MGR00060 
MGR00070 
MGR000N0 
MGRuO»90 
MGR 00100 
MGR 00 1 1 0 
MGR00120 
MGROOl 30 
MGROO Ug 
MGROO) SO 
MuRQOlbO 
MGROOl TO 
MGROO 1 80 
MGROO 190 

mgroo2oo 

MGR00210 
MGR00220 
MGR0U23Q 
MGR00240 
MGR00250 
MGRU0260 
MGR002TO 
MGRUOcNO 
MGR00290 
MGR 00 300 
MGROO JIO 
MGROO 320 
MGR003J0 
MGR 00 340 
MpR00350 
MGROO 360 
MGROO 370 
MGR003NO 
MGHQO J90 
MGR00400 
MGR00410 
MGR004JO 
MGR0O4J0 
MGH00440 
MGP004S0 
MGRU0460 
MGR00470 
MGR00480 
MGR 00490 
MGROOSOO 

mgroosio 
I MGHOOs2o 
MGR00S30 
MGRUOS40 
MGPoOSSO 
MGR00S60 
MGRUOSTO 
MGROOSBO 
MGROOS90 
MGRuObQO 
MGROOblO 
MGR00620 
MGR006 30 
MGR00640 
MGROObSO 
MGRUObbO 
I MGROObTO 
MGR00680 
UGHOO'J'/O 
MGROO TOO 
MGROO MO 
MbRoO 720 
MOP00730 
MUHOO Two 
KjKUO ISO 
mGROO It - 0 
MGROO T To 

m(jm 0 0 mo 

mokuO T90 
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non non non 


ORIGINAL PAGE IS 
OF POOR QUALITY 


tiles MGRIO 


FORTRAN 0 EODL / JOHNSON SPACE CENTEH 


PTIME (MK ) * ( IMIN*1000) ♦IFIX(SEC*10.Q> 
IFIMK.EG.IIGO TO A00 
KCHK*PTIME (MK) -PTIME IHK-l) 
IF(KCHK.GT.30)PTIME(MK)=PTIME(MK)-AQO 
lF(KCHK.LT.O)PTIME(MK)sPTlME (MK) >60000 

A00 continue 

M IS the counter for input file 10 

AO l M*1 

THIS SECTION LOOKS FOR THE BOUNDARIES IN THE OaTA 
DO 3 K*1.NUMB 

IF(MFMRX(K).LT.START(M) )GO TO 3 


L“K ♦ 1 
END»S ' 


• START ( M ) . X W 1 D£ <M) 

DO A N*1 .20 

IF (MFHRx ( L > .GT.ENDIGO TO AA 
A L=L*1 

44 ip!lN&E=MFMRA (LI -MFHRX IK) 

THIS SECTION SETS UP THE FIELD DIMENSIONS 

NPT S=L“K* ] 

ISTART=1 
BTWN*0.0 


THIS SECTION IS FOR LATITUDINAL FLIGHT LINES 

EDGE* «Xw IDE (M)-ARANGE)/2.0 
DO 5 N=1.NPTS 
NUM*L-N> 1 

X(N)=IFIX<ABS(Xn-(EDGE»BTWN) ) ) 

Y(N)*IFIX(YwIDE<M)-MFMRY(NUM) ) 

KTIME(N)=ITIME(NUM) 

DO 6 NN = 1 t IB 

6 MFMR(N.NN)*IMFHR(NUM.NN) 

ktr=l-n 

BTWN=MFMRX (L)-MFMRX(KTR) 

5 CONTINUE 

CALL PRT5(NPTS.ISTART.APHT,ATAT.KT!mE.PTImE.PTEMP.TAT) 

GO TO 55 

THIS SECTION IS FOR LONGITUDINAL FLIGHT LINES 

AS EDGE*(YWIDE(M)-AHANGE)/2.0 
DO 7 N= 1 .NP TS 
NUM=L-N* 1 

X ( N) = I F I X t Xa-MFMRT I NUM > ) 

T (N) =IFIX (ABSIYW- (EDGE*9TWN> ) ) 
ktimein)=itime(num) 

00 B NNsl.lB 

8 MFMR(N.NN) = I MFMR ( NUM . NN ) 

ktr=l-n 

BTWNsMFMPX (L) -MFMRX (KTR) 

7 CONTINUE 

CALL PRT5 (NPTS. I ST ART. APBT.ATAT.KTIME.PT I ME.PTEMP, TAT) 

PREPARE TO WPITE THE OUTPUT FILES 
55 !UNIT=10>m 

•RITE OU. THE FIELDS COMPUTED IN THIS PROGRAM 

IF(ICP.GE.7)G0 TO 333 
00 22 KKsl.NPTS 

MMFMR ( KK ) =Mf MR(KK.l) 

22 WRITE! I UNIT, H 1) IC1. IC2. IC3. IFLLMM! . IL iNt , I RUN . * ( K* ) .T (KK) . 

6 MFMRIKK.I ) .MFMR(KK.2) . (MFMP(KK.MN) ,MN*7,ld) , APR T (KK) .ATAT(KK) 
1 1 1 FORMAT (31 1 .IX. I2.1A.I1,I1. IX. I*. lx. I a.7I1X,1<..I3).1X. <0-6. 2) 

GO TO 28 


MGRUOBOO 

HGR00B10 

MGR00B20 

MGR00830 

MGRUOBAO 

MGRU0B50 

MGR00860 

MGR00d70 

MGR00880 

MGR00890 

MGR00900 

MGR00910 

MGR00920 

MGR00930 

MGR009A0 

MGR00950 

MGR00960 

MGR00970 

MGR00V80 

MGR00990 


MGPO 

MGRO 

MGRO 


MGRO 10 30 


MGRO 

MGRO 

MGRO 

MGRO 


x*=xw 

DE <M) 

MGHO 

y*»yw 

DE IM) 

MGRO 

IF ( |L 

NE.EU.1)XW*0.0 

MGRO 

IF C IL 

NE.EO.A) XW*0.0 

• MGRO 

IF ( IL 

NE.GE.6)YW*0.0 

MGRO 

IF ( IL 

NE.GT.AIG0 TO A5 

MGRO 


MGR01 . 

MGR on 

MGRO 11 
MGRO 1 1 
MGRO 1 j 
MGRO tl 
MGR01200 
MGRO 1210 


OAO 

050 

060 

070 

080 

090 

00 

18 

30 

AO 

50 

60 

70 

80 

90 


MGRO 


MGRO 1230 


MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 

MGRO 


22 0 


2A0 

250 

260 

270 

280 

290 

300 

310 

320 

330 

3A9 

350 

360 

370 

380 


MGRO 1 390 
MGR01A00 
MGP01A10 
MGRO 1 A20 
MGRO 1 A 30 
MGR01AA0 
MGH01A50 
MGR01A60 
MGRO 1a70 
MGRO 1 ABO 
MGRO 1 A90 
MGRO 1500 
MgRO 1510 
mgRO 1 520 
MGRO 1530 
MGR 0 1 SmO 
MGR 01550 
MGR 01560 
MGR01570 
MGR 01580 
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> 1)0 


f/ioe 


QnJGlMAL 

GF FOOH ‘ 


Y 


riLE : MGRID 


FORTRAN 0 EOOL / JOHNSON SPACE CENTER 


J3J CONTINUE 

DO 26 KK = 1 .NPTS 
MMFMR (KK) =mFMH ( KK .3) 

26 WRITE! IiiNIT. 112) I Cl . IC2. IC3. 1 FLO ( M) . IL INE . 1 RUN. X (KK ) . Y (KK ) . 
1 ( MF MR < KK .NM ) ,NM=3,b) .APRT(KK) ,ATAT IHM 

112 FORMAT (31 1. IX. I2.lX.il. 11. IX. 14, lx. 14, 2(lx. 14, IJ), lx. 2F 6. 2) 
26 CONTINUE 


F (M.Eq.nFL'JS) GO TO 33 

M=M* 1 
1LTR*3 

MFLDMFLD(M) 


3 CONTINUE 
33 CONTINUE 
STOP 
ENO 

THIS SUBROUTINE HATCHES THE PRT-5 DATA WITH THE 
CORRESPONDING MFMR DATA 

SUBROUTINE PRT5 (NPTS . I START . APRT . A T AT ,KT I ME ,PT I ME ,PHT . T AT ) 
REAL APR T ( 1 S ) , ATAT ( lbi 
INTEGER PTIME(400) .KTIME(IS) 

REAL PRTI400) .TAT (400) 

KMIN=100000 

kmax=6 

DO V MI=1.NPTS 

IF (KTIME (HI ) .LT . KM IN) KMIN=K TIME (MI ) 

9 IF (KT IME (MI ) ,GT.KMAA)KMAX=KT I ME (Ml ) 

DO 1 J=ISTART.400 
1F(KMIN.GT.PTIME(I))G0 TO 1 


kin. 

m- 1 

GO TO 2 

1 CONTINUE 

2 ISTART=IS 
DO 3 M*LF. 


400 

IF(kmax.gt.ptine(H))go to 3 

LL = M 
GO TO 4 

3 CONTINUE 

4 DO 6 1=1. NPTS 
00 5 N=LF,LL 

IF(KTIME(I; .GT.PTIME(N) IGO TO 5 
M=N-1 

IDIFF=P11hE(N)-PTImE(M) 
K01FF=KTlMt(I)-PTIMt (M) 
DIFFl=KDI r F*1.0 
DIFF2=1DI*F*1.0 
RATI0=01FF1/DIFF2 
DELPRT=PRT (N) -PRT (Ml 
DELTAT=TAT !N) -TAT (M) 

APRT ( I ) =OELPRT»WAT I0*PRT (M) 
ATAT(I)=OELTAT*RATIO*TAT (M) 

GO TO 6 

5 CONTINUE 
b CONTINUE 

RFTURN 
END . 


MGRO 1590 
MGR 01 600 
MGRU1610 
MGhU I o20 
MGR (>1630 
MGR01640 
MGR01680 
MGRU 1660 
MGRO 1670 
MGRU 1680 
MGRO 1690 
MGRU 1 700 
MGHU 1710 
MGRO 1 720 
MGRO 1730 
MGRO 1 740 
MGRO 1 7 SO 
MGRO 1760 
MGRO 1 770 
MGRO 1780 
MGR01790 
MGR01800 
MGRU 1810 
MGRO 1820 
MGRO 1830 
MGRO 1840 
MGH01850 
MGRO 1 860 
MGRO 1870 
MGRO 1 880 
MGRO 1890 
MGRO 1900 
MGRO 1910 
MGRO 1920 
MGRO 1930 
MGRO 1 940 
MGR019S0 
MGRO 1960 
MGRO 1970 
MGRO 1 980 
MGR01990 
MGR02000 
MGR02010 
MGR02020 
HGR02030 
MGR02040 
MGR02050 
MGR02060 
MGR020 70 
MGR02080 
MGR02090 
MGR02100 
MGR021 10 
HGR02120 
MGR02130 
MGR02140 
MGR02150 
MGR02160 
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